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ABSTRACT 



The late-stage evolution of the most massive stars such as r/Carinae is controlled 
by the effects of mass loss, which may be dominated by poorly understood eruptive 
mass ejections. Understanding this population is challenging because no true analogs 
of r\ Car have been clearly identified in the Milky Way or other galaxies. We utilize 
Spitzer IRAC images of 7 nearby (< 4Mpc) galaxies to search for such analogs. We 
find 34 candidates with a flat or rising mid-IR spectral energy distributions towards 
longer mid-infrared wavelengths that emit > 10 5 L in the IRAC bands (3.6 to 8.0 |im) 
and are not known to be background sources. Based on our estimates for the expected 
number of background sources, we expect that follow-up observations will show that 
most of these candidates are not dust enshrouded massive stars, with an expectation of 
only 6±6 surviving candidates. Since we would detect true analogs of r\ Car for roughly 
200 years post-eruption, this implies that the rate of eruptions like r\ Car is less than the 
ccSN rate. It is possible, however, that every M> 40 M star undergoes such eruptions 
given our initial results. In Paper II we will characterize the candidates through further 
analysis and follow-up observations, and there is no barrier to increasing the galaxy 
sample by an order of magnitude. 

Subject headings: stars: evolution, mass-loss, winds, outflows — stars: individual: Eta 
Carinae — galaxies: individual (M33, M81, NGC247, NGC300, NGC2403, NGC6822, 
NGC7793) 



Introduction 



Despite being very rare, massive stars such as luminous blue variable (LBVs), red super giants 
(RSGs), and Wolf-Raye t stars (WRs) play a pivotal role in enriching the interstellar medium (ISM) 
through mass loss (e.g.. lMaederlll98ll ). Understanding the evolution of thes e massive (M>, 30 M ) 
stars is challenging even when mass loss is restricted to continuous winds (e.g.. lFullerton et al.l l2006). 
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but poorly understood im pulsive mass ejections are probably an equally important, if not dom i- 
nant mass loss mechanism ( Humphreys Sz Davidson 1984 ; Smith Sz Owocki 2006 ; Kochaneklboilb ). 
Mass loss also determines the structure of the star at death and hence the observed properties of 
the final core-collapse supernova (ccSN). In addition, there is also evidence that s ome supernova 



(SN) progenitors underg o major mass ejection events shortly before exploding (e.g.. iGal-Yam et al 



2007 



Smith et al.l 120081 ) . further altering the properties o f the explosion and implying a connec 



2007 



Smith k McCrav 



tion between some eruptive mass- loss events an d death (IGal-Yam et al 
20071 ; iKochanekl l2011al ; IChevalier &: Irwinl |2012| ). In two cases, eruptions were observed shortly 



before the ccSN: Type lb SN2006ic was spatially coincident with a bright optical transient that 
occurred in 2004 (jPastorello et al.l 120071 ) . and SN2009ip under went a series of outbursts in 2009 



2010 , and 2011 before proba bly exploding as a Typelln SN (see iMauerhan et al 



2012: 



Pastorello et al.ll2012l ). These processes are likely metallicity dependent (jMevnet et al 



2012; 



Prieto ct al 



1994 



Heger et al.ll2003l ) and there is good evidence that the SNe requiring a dense circumstell ar medium 



(the hyperluminous T ype Iln) predominantly occur in lower metallicity galaxies (e.g., IStoll et al 



2011 



Neill et al 



2011 



)■ 



Traditional studie s of these massive stars search for them optic ally and then characterize them 



spectroscopically (e.g.. iBonanos et al.ll2009i . 



2010 



Clark et al.ll2012l ). This approach is not ideal for 



probing the episodes of major mass- loss because of dust formation in the ejecta. Dense winds tend 
to form dust, although for h ot stars the wind must be dense enough to form a pseudo-photosphere 
in th e wind ( Davidson! 119871 ) that shields the dust formation region from the UV emission of the 
star (IKochanekl l2011bl ). Th e star will then be heavily ob scured by dust for an extended period 
after the eruption (see, e.g., Humphreys &, Davidson! Il994l ). The great eruption of 7/ Car between 
1840 and 1860 is the most famous case of a st ellar outburst, ejecting ~ 10M Q material before 
reappearing as a hot star in the 1950s (see, e.g.. Humphreys et al.lll997l ). The ejecta are now seen 
as a dusty nebula around the star absorbing and then reradiating ~ 90% of the light in the mid-IR. 
This means that dusty ejecta are a powerful and long-lived signature of eruption. The emission 
from these dusty envelopes peaks in the mid-IR with a c haracteristic red c olor and a rising or flat 
spectral energy distribution (SED) in the Spitzer IRAC (jFazio et al.ll2004l ) bands. 



In the Galaxy, s t ars with resolved shells o f dust emission primarily at 24 ^m are easily found 



(jWachter et al.l l2010l : iGvaramadze et ajj l2010l ) . The advantage of the 24 [im band is that it can 
be used to identify dusty ejecta up to 10 3 — 10 4 years after its formation. A minority of these 

10 55 LfTi) with massive (~ 0.1 — 10 Mr) shells (see sum- 



2001 



Humphrevs et al. 


1999: 


Smith & Owocki 


2006; 


Smith 


(Voors et al. 


20001. the Pistol Star ( 


Figer et al. 


1999 





GjU29+0.46 (jffiggs et al1l994h . Wray 17-96 (fcgan et alhooj ). and IRAS 18576+0341 (jlJeta et al 



These systems are far older (> 10 3 years) than r\ Car, which makes it difficult to use the ejecta 
to probe the rate or mechanism of mass-loss. Still, the abundance of Galactic she lls implies that th e 
rate of r] Car- like eruptions is on the order of a modest fraction of the ccSN rate ()Kochanekll2011bl ). 
Their emission peaks in the shorter IRAC bands when they are relatively young (~ 10 — 100 years) 
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because, as t he ejected material expa nds, the dust becomes cooler and the emission shifts to longer 
wavelengths ([Kochanek et al.l l2012al ). It is difficult to quantify searches for such objects in our 
Galaxy as it is difficult to determine the distances to the sources and the survey volume because 
we have to look through the crowded and dusty disk of the Galaxy. Surveys of nearby galaxies are 
both better defined and build larger samples of younger systems whose evolution can be studied to 
better understand the mechanism. 

With Spitzer it is difficult to use the 24 /jm observations that have proved so successful in the 
Galaxy because of the poor angular resolution. However, w e hav e shown tha t such surveys can 
be done with IRAC (3.6-8.0 /im). In iThompson et al.l (|2009l ) and iKhan et al.l (|2010i ). we charac- 
terized the extreme AGB sta r populations that appear to be the progenitors of the SN2008S-like 



transients (Prieto et al 



2008 



Prietol 120081 ) using archival IRAC images of nearby galaxies. These 
studies empirically confirmed that these ~ 10 4,5 L Q dusty stars are rare but are also relatively 
easy to identify in IRAC images despite the mode st angular resolutio n. Next, we examined all the 
other bright, red mid-IR sources in M33, and in IKhan et al.l (|201ll ) we discovered Object X, the 
brightest mid-IR star in M33. Object X is a Ly ~ 5 x 10 5 L Q , M> 30M Q evolved star obscured by 
dust formed during mass loss events over the last ~ 1 century Its properties are similar to those of 



the G alactic OH/IR star IRC+10420 (jHumphrevs et al.lll9971 ; lBlocker et alJll999 : Humphreys et al 



2002), which has a complex dusty circumstellar structure resulting from episodic, low- velocity mass 
ejections. We proposed that Object X may emerge fro m its current ultra-short e volutionary phase 



as a hotter post-RSG star analogous to M 33 Var A (jHubble fc Sandage 
200fil ). 



1953 



Humphreys et al 



While Object X is intriguing, it likely underwent a period of enhanced, but relatively steady, 
mass loss from the parent star rather th an the short transient episode of mass loss usually associated 
with so called "supernova impostors" (|Van Dvk et al.ll2000l ). It is also an order of magnitude less 
luminous and several times less massive than 77 C ar, one of the most luminous ( Lhni. ~ 5 x 10 6 Lq) 
and massive (M ~ 100 — I5OM0) stars known (jHumphrevs fc Davidson! I1994T). No true analog 



of n Car in mass, lumi nosity, energetics, mass lost and age has been found (see lSmith et al.l 12011 



Kochanek et al.ll2012al ). Quantifying the population of 77 Car analogs, or their rarity, in the local 
universe can allow us to investigate the rate of giant eruptions of the most massive stars. It may 
also help us answer open questions about the evolution of massive stars such as: (1) the frequency 
of major mass ejection events, (2) the number of events per star, (3) whether the frequency depends 
on the metallicity or other stellar properties, and (4) whether there is really any relation between 



mass ejections and the so called "supernova impostors" (see ISmith et al.l l201ll ; Kochanek et al 
2012ah . 



Here we carry out a pilot study of 7 nearby galaxies within (< 4 Mpc) to search for analogs of 
77 Car. We concentrate on galaxies with recent star formation, as only these will have large numbers 
of the short-lived, very massive stars that we want to study, but we also include one small, low-mass 
galaxy (NGC6822) as a test case (see TableU]). Section[2] describes our methodology for identifying 
potential r\ Car analogs in nearby galaxies using archival Spitzer data and sources of contamination. 
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Table 1: Properties of Targeted Galaxies 



Galaxy 


Distance 


M B 


E(B - V) 


log 1Q L{Ha) 


SFR (Ha) 


Known ccSN 




(Mpc) 






(erg/s) 


(M /yr) 


(< 20 years) 


NGC 6822 


0.46 


-14.9 


0.24 


39.1 


0.01 




M33 


0.96 


-18.8 


0.04 


40.6 


0.33 




NGC 300 


1.9 


-17.7 


0.01 


40.1 


0.11 




NGC 2403 


3.1 


-18.7 


0.04 


40.8 


0.44 


SN2004dj (IIP) 


M81 


3.6 


-20.1 


0.08 


40.8 


0.46 


SN 1993J (lib) 


NGC 247 


3.6 


-18.2 


0.02 


40.3 


0.17 




NGC 7793 


4.1 


-18.5 


0.02 


40.6 


0.33 


SN2008bk (IIP) 



Section[3] discusses the nature of the candidates, although a detailed study is deferred to Paper II. 
Section^] shows how our search method allows us to quantify the selection criteria and to set an 
interesting limit on the rate of events similar to the Great Eruption of r] Car in the local universe 
even before we have completed Paper II. Finally, in Section[S]we outline the future of our approach. 



2. A Search for r/ Car Analogs 

In this section, we present the methodology of our search for rj Car analogs. First we discuss 
our data sources and the properties of the targeted galaxies. Next we describe the photometry and 
candidate selection procedures. Then we consider contamination due to non-stellar sources. Finally 
we consider if an r\ Car analog in a nearby galaxy could be hidden in a compact stellar cluster. 



2.1. Targeted Galaxies 



There are a number of sources of archival Spitzer data for nearby galaxies. In the Local 



Group, SAGE (Surveyin g the Agents of a Galaxy's E volution) and SAGE- SMC (IMeixner et al 



which allowed for mid-IR variability studies of M 33 stars (McQuinn et al 



20061 : iGordon et al.ll2007h surveyed the L MC and SMC Barmbv et al.1 (|2006j) surveyed M31 with 
first results for massive stars discussed by I Mould et al.l (120081) . M 33 was observed at several epochs 



2007 



Thompson et al 



2009). The Spitzer Infrared Nearby Galaxies Survey (SINGS. iKennicutt et al.ll2003l ) made a com- 
prehensive mid-IR imaging and spectros copic survey of 7 5 galaxies, many of them within 10 Mpc. 
The Local Volume Legacy Survey (LVL, iDale et al.ll2009i ) surveyed a total of 256 nearby galaxies, 
including all known galaxies inside a sub- volume bounded by 3.5 Mpc and an unbiased sample of 
S-Irr galaxies within a larger, and more r epresentative, 11 M pc sphere. The ongoing Spitzer Sur- 
vey of Stellar Structure in Galaxies (5 4 G, ISheth et alJl2008l ) is collecting data for ~ 2300 galaxies 
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within 40 Mpc using the warm Spitzer (3.6 and 4.5/xm) bands. 



For our paper, we selected 7 galaxies spanning a range of mass, morphology, distance, and star 
formation history. Since this is a pilot study, the sample is deliberately eclectic rather than focused 
on a sample maximizing the star formation rate per galaxy. Ultimately we would like to examine 
all nearby galaxies rather than just a few. TableU] summarizes the pr operties of the targete d 
galaxies. The absolute magnitude Mg and Ha luminosity L(Ha) are from iKennicutt et al.l ( 2008). 
and L(Ha) is converted to star formation rate (SFR) following Equatio n 2 of lKennicutt] (|l998l ). The 
foreground Galactic extinctions are from lSchlaflv &: Finkbeinerl (|201ll ). The targete d gala xies have 
an integrated SFR of ~ 2 M© year -1 . For the assumed Salpeter IMF of lKennicutt! ( 19981 ). we can 
convert this to the massive (M > 8 Mq) star formation rate of ~ 0.014 year -1 . The observed ccSN 
rate over that 20 years is ~ 0.15 year -1 (0.05 < Rsn < 0.35 year -1 , at 90% confidence). Since the 
ccSN rate should agree with the massive-star formation rate, this is a significant discrepancy for 
which we have no ob vious explanation, and such mismatches are also found in other contexts (e.g., 
Horiuchi et al. 2011). 



The nearest of our targeted galaxies, NG C 6822 (D ~ 0.46 Mpc. lGieren et al.l l2006). is a barred 
irregular galaxy (jde Vaucouleurs et al.lll99ll ). We included this small galaxy in our sample as an 
interesting nearby t est case for examini ng large numbers of smaller, lo wer metallicity systems . 



M3 3 (D~ 0.96 Mpc. LBonanos et al.ll2006l ) was previously studied by both iThompson et al.1 1)20091 ) 



and iKhan et al.1 (|2010l ) to search for dusty stars that are much redder but less luminous than 
the stars we are se arching for in thi s paper. NGC 300 (D ~ 1.9 Mpc, iGieren et al.l 120051 ) and 



M81 (D ~ 3.6Mpc, ICerke et allhoilh were also studied bv khan et all (20ld). NGC 2403 (D 



3.1 Mpc, ISaha et al l hood ) contai ns two sources sometimes classified as SN impostors, SN 1954J 



and SN 2002kg (see the review bv lvan Dvkll2005l ). but any star associated with SN1954J must be 
relatively low mass (~ 20 M© rather than > 50 M©) and shows no strong evidence for mid-IR emis 



Kochanek et al 



sion, w hile SN 2002kg is simply a luminous variable star with little mass lo ss (see 
201231 ). Unlike the other large galaxies w e studied, NGC 2 47 (D ~ 3.6 Mpc. iMadore et al.ll2009h is 
highly inclined. NGC 7793 (D ~ 4.1 Mpc. ItuIIv et allhood ) is the most distant galaxy studied. 



For M 33 , we used the six co- ad ded epochs of IRAC da ta from iMcQuinn et al.1 (|2007l ) that 
were used bv IThompson et al.l (|2009l ) and IKhan et al.l (|2010l ). and the MIPS data retrieved from 
the Spitzer Herit age Archive. For NGC 300 and NGC 247, we used the data from the LVL survey 
(|Dale et alJ boosh. For NGC 6822. N GC 2403. M81, and NGC 7793, we used the data from the 
SINGS survey (jKennicutt et alj|2003l ) . We utilize the full mosaics available for each galaxy. Table[2] 
shows the different pixel scales of the images retrieved from the Spitzer archive for M 33, and those 
provided by the SINGS and LVL surveys for the other six galaxies. 
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2.2. Candidate Selection 

The SED of a hot dust-obscured star will generally have two peaks — a dust obscured optical 
peak, which could be absent altogether given enough absorption, and a mid-IR peak whose location 
in wavelength depends on the radius of the dust shell around the star. In the IRAC bands, the 
SED will be flat or rising towards longer wave lengths. For example , 77 Car has a steeply rising 



SED towards longer mid-IR wavelengths (e.g., iRobinson et al.1 Il973l ) and the luminosity of the 



star exceeds 10 5 Lq in each IRAC band (see Figure[T]). At extra-Galactic distances, an 7/ Car 
analog would appear as a bright, red point source in IRAC images with a relatively fainter optical 
counterpart due to the self-obscuration. 



We used the Daophot/ Allstar PSF-fitting and photometry package (|Stetsonlll992l ) to identify 
point sources in all four IRAC bands and then performed photometry at the source location using 
both aperture and PSF photometry. We used the IRAF0 ApPhot/Phot tool for the aperture 
photometry. The aperture fluxes were transformed to Vega-calibrated magnitudes following the 
procedures described in the Spitzer Data Analysis Cookbook and aperture corrections of 1.213, 
1.234, 1.379, and 1.584 for the four IRAC bands. The choice of extraction aperture aperture (R ap ) 
as well as the inner (Ri n ) and outer (R ou t) radii of the local background annulus are reported in 
TableO We estimate the local background using a 2a outlier rejection procedure in order to exclude 
sources located in the local sky annulus, and correct for the excluded pixels assuming a Gaussian 
background distribution. Using a background annulus immediately next to the signal aperture 
minimizes the effects of background variation in the crowded fields of the galaxies. We used the 
Daophot /Allstar package for PSF photometry. The PSF photometry fluxes were transformed to 
Vega-calibrated magnitudes by applying zero point offsets determined from the difference between 
the calibrated aperture magnitudes and the initial PSF magnitude estimates of the bright stars in 
each galaxy. 

For the 3.6 and 4.5/jm bands, after verifying consistency with the aperture magnitudes, we 
only use the Vega-calibrated PSF magnitudes. For 5.8 /im, we switch to aperture magnitudes when 
Allstar fails to fit the PSF to a point source at the location identified by Daophot due to the 
decreasing resolution. PSF photometry performs very poorly at 8.0 /jm, leading to both inaccurate 
photometry and many false sources because Daophot frequently splits up extended regions of PAH 
emission into spurious point sources. Thus, at 8.0 /im we only use aperture photometry at positions 
determined for sources identified in the other three bands. We do not use this band for building 
our initial source list. 

We define our initial source list as all point sources that have \L\ > 10 4 Lq in any one of the 
3.6, 4.5, and 5.8 bands, excluding regions near saturated stars and, in the case of M81, the 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation. 

2 http://irsa.ipac. caltech.edu/data/SPITZER/docs/dataanalysistools/ 



-7- 



Table 2: Aperture Definitions 



Band 


Pixel Scale 


Rap 


Rin 


Rout 


(pm) 


(Archive) 


(Survey) 








3.6-8.0 


1"2 


0"75 


2"A 


274 


7"2 


24 


2"45 


1*5 


375 


670 


8"0 


70 


4"0 


4*5 


1670 


1870 


3970 


160 


8"0 


9"0 


1670 


6470 


12870 



high surface brightness core of the galaxy. We identify sources in each of these three bands, and 
cross-match the catalogs using a 1 pixel matching radius. We then adopt the position determined at 
the shortest wavelength (highest resolution) with a > 3a detection, and we use this position for the 
8.0 fim aperture photometry. We fit the mid-IR SED of each object as a power law in wavelength 

log 10 (AL A ) = a x log 10 (A) + b (1) 

to determine the slope (a, \L\ oc A a ) and intercept (b). We can crudely relate the slope (a) to a 
dust temperature as 



log 



10 



-4 + (2) 

where Ai and A4 are the shortest and longest band-centers assuming a blackbody spectrum and 
ignoring Planck factors. We define the total mid-IR luminosity (L m in) as the trapezoid rule integral 
of L\ across the band centers 

3 1 

L m iR = ^2 - (Ai+i - Aj) (L\. + L Xl+1 ) , (3) 
i=i 

where Aj = 3.6, 4.5, 5.8, and 8.0 /mi. We also calculate the fraction / of L m jR that is emitted in 
the first three IRAC bands. We define / as the ratio of the energy emitted between 3.6 and 5.8 fim 
(first two terms of the integral), to L m m (all three terms of the integral). The approximate values 
of L m m, a, and / for ryCar are 10 5 ' 65 L , 2.56, and 0.32, and those for Object X are 10 517 L , 
0.22, 0.57. 

We defined candidates as sources with mid-IR luminosity L m m > 10 5 L Q , a mid-IR SED 
slope a > 0, and / > 0.3. Figures [2] and [3] show the distribution of point sources in M81 with 
\L\ > 10 4 L in at least one of the 3.6, 4.5, and 5.8 /xm IRAC bands as a function of L m jR, a, and 
/. The open red triangles in these figures correspond to candidates that are known to be non-stellar 
in nature (see Section [2T3|) . and the solid red triangles represent the surviving candidates. While a 
few hundred sources in M81 are bright enough in the mid-IR to be included in these figures, only 
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a handful of these even remotely resemble rj Car, and not a single one of them is as luminous and 
as red (cold) as ij Car. The other targeted galaxies show similar distributions of sources. These 
distributions illustrate that our selection criteria for identifying potential 77 Car analogs are robust 
and allows for selecting objects that are significantly less luminous in the mid-IR and have much 
warmer circumstellar dust than r\ Car. Table[5]reports the survey area and the number of candidates 
found for each galaxy. 

We used aperture photometry to estimate the MIPS 24, 70, and 160 ^m band luminosities of 
the objects that meet our selection criteria. For point sources that do not have a flux that is > 3a 
above the local sky, we determine the 3a detection limit for each aperture location using the local 
background estimate. Due to the poor spatial resolution of these bands, which forces us to choose 
increasingly large apertures at longer wavelengths (see TableE]), these measurements have limited 
utility. Figured shows the mid-IR SEDs of the candidates we identified in M33 along with normal 
stars in the M 33 image selected from top left region of Figure[5j At 24 /iin, the SEDs of the normal 
stars show the expected slope for the Rayleigh-Jeans tail of their SEDs, followed by an unphysical 
rise at 70 and 160 fim. Essentially, due to the poor resolution, the apertures used for these two 
bands include many objects other than the intended target, and even normal stars appear to have 
rising far-IR SEDs. This means that we can generally use the 24 /um fluxes while the 70 and 160 /xm 
measurements should be treated as upper limits regardless of their origin. Nevertheless, the MIPS 
bands are useful as a qualitative constraint on an object's physical nature (i.e. if it is a galaxy, 
QSO, cluster etc.). 



2.3. Sources of Contamination 



While our selection criteria are designed to identify dust obscured individual stars, there are 
sever al classes of conta minating sources as well. QSOs have red mid-IR SEDs compared to stars 



(e.g., IStern et a" 



SED models in 



2005), as do sta r forming galaxies with strong PAH emission at 8.0 /jm (e.g. the 
Assef et al.l l2010l ). Sources in the galaxies such as dusty star clusters and H II 
regions can also appear as candidates. We used the SIMBAD^I and VizieR0 services to search for 



previ ous classifications and near-IR counterparts from the 2MASS point source catalog (jCutri et al 



2003). We also noted other significant pieces of information, such as a candidate location coinciding 
with known radio and X-ray sources. For example, many of the candidates in M33 are associated 
with known supernova remnants (SNRs). We reject a source as "non-stellar" if it is a galaxy or 
QSO with a measured redshift or if archival images clearly show that it is a galaxy. The rejected 
sources that meet our selection criteria are described in Section l3~Tl 



We estimate the expected number of extragalactic contamination for each galaxy using the 



3 http://simbad. u-strasbg.fr/ 
4 http: //vizier, u-strasbg.fr/ 
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Table 3: Candidate Statistics 





NGC 


M33 


NGC 


NGC 


M81 


NGC 


NGC 




6822 




300 


2403 




247 


7793 


Survey Area A (deg 2 ) 


0.1 


0.73 


0.17 


0.12 


0.17 


0.2 


0.044 


Candidates 





9 


1 


5 


14 


3 


14 


Expected Background £ (deg -2 ) 


0.17 


1.67 


11 


47 


75 


75 


110 


Expected Contamination AT, 





1 


2 


6 


13 


15 


5 


Rejected Candidates 














7 


1 


3 


Remaining Candidates 





9 


1 


5 


7 


2 


11 



SDWFS survey (jAshby et al.ll2009l ) where th e nature of the sources, particularly AGNs, is also well 
understood from the AGES redshift survey (jKochanek et al.ll2012bl ). We transform the apparent 
magnitudes of all sources in a 6 deg 2 region of SDWFS to luminosity using each target galaxy's 
distance modulus, determine how many of them would meet our selection criteria, and correct 
that count for our survey area around each galaxy. Table[3] reports the expected surface density 
of extragalactic contaminants and the number expected given the survey area around each galaxy. 
We expect a total of ~ 41 extragalactic sources to pass our selection criteria across the targeted 
galaxies, as compared to 46 initial candidates. FigureHJ which has the same format as Figure[2j 
illustrates this for M81's distance. In the 6 deg 2 SDWFS area, 449 (~ 75 deg 2 ) sources pass our 
selection criteria, indicating that we should expect ~ 13 background sources meeting our selection 
criteria given our 0.17 deg 2 survey region around M81, as compared to the 14 initial candidates 
selected. Indeed, as we discuss in SectionHl we can already identify 11 of the 46 initial candidates 
as extragalactic. Statistically, this means that only 6 ± 6 are likely associated with the galaxies. 
Also note in FigureH] that none of the contaminating background sources have properties directly 
comparable to r\ Car. The expected numbers of contaminating sources are generally consistent with 
the observed numbers with the exception of NGC 247, which we investigated but appears to surely 
be a statistical fluke. The angular distribution of the candidates relative to the galaxies is also 
strongly suggestive of a dominant contribution for background sources. 

Many of the candidate SEDs show a "dip" from 3.6 \xm. to 4.5 /im before rising again at 5.8 /im 
(see FigureE l). This is a co mmon feature of star cluster SEDs created by strong PAH emission at 



3.6 /im (jWhelan et al 



2011 



L). In total galaxy spectra, this is a weaker effect and the dominant PAH 
emission feature is at 8/um and comes more from the diffuse ISM rather than individual stars or 
clusters. The SEDs of 77 Car and Object X do not show this dip at 4.5 /im. We treat the presence 
of this dip as a qualitative indicator that the source may be a cluster or lie in a cluster. Deep 
Hubble Space Telescope (HST) images of these regions, where available, can help us distinguish 
single bright red stars from clusters of fainter stars that may be merged into a single bright source 
in the lower resolution Spitzer images. Whether these clusters can potentially hide r\ Car analogs 
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is discussed in Section [2~4l 

In Figure[71 we present SEDs of four different types of objects that met our selection criteria 
— a likely dusty star in NGC 2403, a star-cluster in M 33, a QSO behind M 81, and a galaxy behind 
NGC 7793. Although all four objects met our selection criteria, the detailed shapes of their SEDs 
are very different from each other. The likely stellar source, N 2403-3, has a very steeply rising 
mid-IR SED that peaks between 8 //m and 24 /im. While the compact cluster (M 33-8) SED looks 
quite similar to that of 77 Car in the IRAC bands, it continues to rise steeply up to 24 //m, and the 
MIPS 70 /im and 160 /mi upper limits show that it peaked between 24 /mi and 70 /mi. The SED 
of the QSO (M81-4) remains relatively flat from 3.6 /mi to 24 /mi. The MIPS 70 /im and 160 /im 
upper limits for N 7793-2 (a galaxy) are quite stringent, because the source is far from the center 
of the galaxy, and would rule out an 77 Car analog model. For some cases, such as N 7793-2, HST 
images clearly determine the nature of the source (FigurefTOl bottom panel). 



2.4. Star Clusters 



One concern with star clusters as a source of contamination is the possibility of "hiding" a 
luminous dusty star in a dusty star cluster. To explore this problem we estimated what the SED of 
the star cluster containing r\ Car would look like if it were located in one of the targeted galaxies. We 
comb ined the SED of 7 7 Car fr om lHumphrevs Davidson! (|1994l ) with the SED of the Carina nebula 
from ISmith &: Brooks! (|2007l ) to produce an SED of the entire complex (Figure[6]) . The combined 
SED is clearly a multi-component SED, which is not typical of our candidates. Moreover, the 
Carina nebula is roughly ~ 2.5° in extent (ISmith &: Brooks! 120071 ) . which at the distance of M33 
becomes ~ 20'/ and would be easily resolved by IRAC. Even at the distance of NGC 7793 it would 
still subte nd ~ h". and be resol ved. At all these distances it would be very easily resolved by HST 
or JWST ([Gardner et allboOfih . 



In Figures [2] and [3l we show the mid-IR luminosity L m m, SED slope o, and fraction / of 77 Car 
('V')j the Carina nebula excluding 77 Car ("77— "), and the entire complex including 77 Car ("77+"). 
It is apparent from these figures that even if the Carina nebula was not resolved: (1) we would 
select analogs of 77 Car and unresolved dusty stellar complexes hosting such analogs, (2) while it is 
close, we would not select a stellar complex that is identical to the Carina nebula excluding 77 Car, 
and (3) there are no sources with L m m, a, and / comparable to 77 Car in M81. Indeed, this last 
point is true for each galaxy we studied. 

There are, however, far more compact star clusters among the candidates such as M 33-5, M 33- 
8 and M 81-10 (see Section[3j) where HST images are required to recognize their spatial extent. Even 
in these cases it is unlikely we would lose a candidate. First, it would require a "conspiracy" of 
a sort, namely that the SED of the hotter circumstellar dust around the star (with characteristic 
T ~ 400 K and L^oi ~ few x 10 6 Lq) seamlessly merges with the colder SED of the interstellar dust 
(with characteristic T ~ 100 K and Lboi ~ 10 7 Lq) in the cluster. Typically we find that this leads 



Table 4: Rejected Candidates 



ID 


KA 


Dec 


blope 


lo SlO L 


£ 

J 


Spectral Energy Distribution [lo 


5 10 (AL A /L Q )] 








(a) 




[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


[70] 


[160] 


M81-1 


149.21474 


69.12843 


0.56 


5.50 


0.55 


5.46 


5.55 


5.61 


5.66 


5.46 


<5.28 


<5.81 


M81-2 


148.91331 


69.29649 


0.22 


5.06 


0.54 


5.21 


5.15 


4.95 


5.33 


5.13 


<5.53 


<5.97 


M81-3 


149.25626 


68.91674 


0.66 


5.43 


0.53 


5.37 


5.45 


5.57 


5.59 


5.58 


<5.50 


<5.60 


M81-4 


149.15222 


69.00780 


0.61 


5.08 


0.53 


5.08 


5.09 


5.19 


5.28 


5.30 


<5.63 


<6.58 


M81-8 


149.17365 


68.80576 


1.14 


5.42 


0.47 


5.30 


5.32 


5.57 


5.66 


5.60 


<5.51 


<5.60 


M81-9 


148.85034 


69.24746 


1.76 


5.45 


0.40 


5.18 


5.38 


5.45 


5.83 


5.52 


<6.46 


<6.27 


M 81-13 


148.99657 


69.26839 


1.74 


5.03 


0.31 


4.90 


4.85 


4.78 


5.54 


5.10 


<6.14 


<5.99 


N 247-2 


11.90806 


-20.51950 


1.35 


5.10 


0.36 


5.13 


4.85 


5.01 


5.55 


5.46 


<6.89 




N 7793-2 


359.47302 


-32.47820 


0.62 


5.37 


0.51 


5.39 


5.33 


5.49 


5.57 


5.84 


<5.99 


<6.00 


N 7793-5 


359.34467 


-32.62253 


0.76 


5.15 


0.42 


5.33 


4.96 


5.05 


5.54 


5.24 


<6.13 


<6.47 


N 7793-7 


359.46878 


-32.63801 


1.17 


5.04 


0.37 


5.09 


4.87 


4.87 


5.49 


5.59 


<6.08 


<6.80 



to SEDs with "bumps" which we do not observe. 



Possibly more constraining is the requirement that for a compact cluster to hide an rj Car 
analog it must still contain large amounts of interstellar gas and dust several million years after 
the cluster formed to allow for the time that even the most massive stars require to evolve away 
from the main sequence. However, a cluster sufficiently luminous to hide an rj Car analog must host 
many luminous stars with strong UV radiation fields and winds, which will likely clear the cluster 
of gas and dust needed to produce strong mid-IR emission. For example, 30 Dor, which harbor s 
stars possibly as massive as 300 M Q and is about 1.5 x 10 6 years old (e.g.. ICrowther et al.l l2010i ) , 
is a weak source of 8/^m emission (see, e.g., Figure 1 of lZhang fe StanekH2012l ). 



3. Inventory of Candidates 



In this section we discuss the initial results of our search for 77 Car analogs. A total of 46 sources 
passed our basic mid-IR selection criteria {L m m > 10 5 L , a > 0, / > 0.3). First we discuss the 
eleven candidates that can be rejected as known non-stellar sources. Their properties are reported 
in Tabled! Then we list the remaining 35 candidates, including Object X, in Table[5j Table[6] 
presents the near-IR phot ometry for the 9 sources with counterparts in the 2MASS point-source 
catalog (jCutri et al.l l2003). 



-12- 



3.1. Rejected Candidates 



Of the 11 rejected candidates, six are AGNs or galaxies with a redshift measurement . Four have 
been photometrically classified as galaxies by the SDSS survey (jAbazaiian et al.ll2009l ) and visual 
inspections of the SDSS images find diffuse extended sources consistent with this classification. One 
of these sources is also 0'.'4 away from a radio and X-ray source and is likely a low redshift AGN. 
One candidate is a well-resolved galaxy in HST images. The SEDs of the rejected candidates are 
shown in Figured] and their luminosities, SED slope, and / are listed in Tabled! In detail, we find 
that: 



M81-1 is an AGN. It is 0"362 from a quasar at z = 0.605 (jBichards et al.l 12003 ). 



M 81-2 is an AGN. It lies Of! A from a radio and X-ray source (|Fleschll2010i ) and is classified as 
a galaxy by SDSS (jAbazajian et al.ll2009l ). Visual inspection of the SDSS image also shows 
an extended source consistent with this classification. 



M81-3 is an AGN. It is 0"38 from a quasar at z = 1.29683 ((Schneider et al.ll2010h . 



M81-4 is an AGN. It is Of! 65 from a quasar at z = 1.97519 ((Schneider et al.l l20ld ). 



M81-8, M81-9, and M 81-13 are classified as galaxies by SDSS (jAbazajian et ajJ|20Q9j) and 
visual inspections of the SDSS images finds extended sources consistent with these classifica- 
tions. 



N 247-2 is a galaxy at z = 0.02089 (| Jones et al J 12003 ). 

N 7793-2 is unambiguously a galaxy in HST images (Figure[T0| bottom panel). 



N 7793-5 is a galaxy at z = 0.0614 (I Jones et alJl2009l ). 



N 7793-7 is an AGN. It is 1'.'5 from a QSO at z = 0.071 (jVeron-Gettv fc Veronl 12010 1. 



While the astrometric matches are sometimes imperfect, we are dealing with objects with 
such low surface densities that a mismatch is e xtraordinarily unlikely. Essentially, this search 
recapitulates aspects of the iKozlowski et al.l (120101 ) search for quasars behi nd the Mage l lanic clouds 
as red mid-IR sources following the extragalactic mid-IR search criteria of lStern et al.l (|2005l ). The 
red mid-IR colors created by the power-law SEDs of quasars mimic aspects of the red SEDs of 
dusty stars. 



3.2. Remaining Candidates 



The mid-IR luminosities, slopes, and the fractions / of the 35 remaining candidates, including 
Object X, are presented in Table[5j Of these, 9 also have 2MASS photometry and their near- 
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Table 5: Remaining Candidates 



ID RA Dec Slope i ogl0 %p / Spectral Energy Distribution [log 10 (AL.\/ L©)] 













( 


a) 










[3.6] 


[4.5] 


[5 


.8] 


[s 


.0] 


[24] 


[70] 


[160] 


M33-l a ' 6 


23 


35015 


30 


42626 





24 


5 


11 





57 


5 


14 


5 


17 


5 


25 


5 


21 


5 


13 


<5 


78 


<6.09 


M 33-2 a 


23 


39209 


30 


69071 


1 


06 


5 


10 





46 


5 


06 


4 


91 


5 


27 


5 


34 


6 


10 


<6 


90 


<6.76 


M 33-3 


23 


43939 


30 


61357 


1 


98 


5 


10 





39 


4 


87 


4 


69 


5 


31 


5 


42 


5 


51 


<6 


15 


<6.13 


M 33-4 


23 


55650 


30 


56175 


2 


17 


5 


01 





36 


4 


77 


4 


52 


5 


21 


5 


36 


5 


64 


<6 


44 


<6.58 


M 33-5 


23 


31891 


30 


88054 


2 


32 


5 


57 





36 


5 


19 


5 


33 


5 


67 


5 


90 


6 


86 


<6 


90 


<6.54 


M 33-6 


23 


36988 


30 


67363 


1 


79 


5 


05 





36 


4 


90 


4 


81 


5 


08 


5 


47 


6 


20 


<6 


57 


<6.33 


M 33-7 


23 


39793 


30 


65805 


2 


33 


5 


08 





36 


4 


76 


4 


63 


5 


28 


5 


43 


5 


73 


<6 


33 


<5.99 


M 33-8 a 


23 


50089 


30 


67987 


2 


40 


5 


56 





35 


5 


22 


5 


13 


5 


72 


5 


94 


6 


46 


<6 


72 


<6.50 


M 33-9 


23 


37096 


30 


67276 


2 


34 


5 


02 





32 


4 


73 


4 


61 


5 


10 


5 


45 


6 


12 


<6 


56 


<6.34 


N 300-1 


13 


71123 


-37 


67159 





96 


5 


27 





50 


5 


17 


5 


27 


5 


37 


5 


50 


5 


28 


<6 


23 


<6.30 


N2403-P 


114 


32964 


65 


59473 


1 


37 


5 


22 





38 


5 


20 


5 


03 


5 


20 


5 


63 


5 


96 


<6 


96 


<7.37 


N 2403-2 


114 


20582 


65 


60922 


2 


00 


5 


11 





36 


4 


93 


4 


68 


5 


25 


5 


49 


5 


31 


<7 


03 


<7.57 


N 2403-3 


114 


09702 


65 


61411 


2 


56 


5 


16 





33 


4 


75 


4 


83 


5 


24 


5 


59 


5 


53 


<6 


77 


<7.33 


N 2403-4 a 


114 


22210 


65 


59257 


2 


20 


5 


17 





32 


4 


98 


4 


71 


5 


24 


5 


61 


5 


24 


<6 


71 


<7.45 


N 2403-5 


114 


22632 


65 


59669 


2 


14 


5 


20 





30 


4 


98 


4 


89 


5 


15 


5 


69 


5 


44 


<6 


90 


<7.50 


M81-5 


148 


75421 


69 


12405 





61 


5 


01 





52 


5 


01 


4 


99 


5 


13 


5 


20 


5 


22 


<6 


06 


<7.12 


M81-6 


148 


83128 


68 


95947 





89 


5 


01 





52 


4 


88 


5 


02 


5 


15 


5 


19 


5 


08 


<5 


50 


<6.56 


M81-7 


148 


72035 


69 


14713 





73 


5 


05 





51 


5 


03 


5 


07 


5 


11 


5 


29 


5 


31 


<6 


74 


<7.07 


M 81-10 


148 


97075 


68 


98440 


1 


72 


5 


57 





37 


5 


39 


5 


43 


5 


54 


5 


99 


6 


25 


<7 


20 


<7.43 


M 81-11 


149 


00545 


68 


98351 


1 


71 


5 


10 





37 


4 


95 


4 


89 


5 


12 


5 


51 


5 


20 


<6 


70 


<7.46 


M 81-12 


149 


01479 


68 


98553 


2 


46 


5 


05 





32 


4 


71 


4 


65 


5 


15 


5 


47 


5 


45 


<6 


70 


<7.46 


M 81-14 


148 


66461 


69 


08003 


2 


63 


5 


08 





30 


4 


74 


4 


57 


5 


19 


5 


53 


5 


05 


<6 


92 


<7.28 


N 247-1 


11 


51449 


-20 


72367 


1 


12 


5 


14 





50 


4 


97 


5 


15 


5 


25 


5 


38 


5 


28 






<5.36 


N247-3 a 


11 


91676 


-20 


90363 


1 


45 


5 


14 





31 


5 


07 


5 


02 


4 


64 


5 


68 


5 


32 


<5 


95 


<6.33 


N 7793-1 


359 


39191 


-32 


54715 





31 


5 


13 





55 


5 


20 


5 


28 


4 


95 


5 


40 


5 


21 


<5 


45 


<6.47 


N 7793-3 


359 


43268 


-32 


60958 





77 


5 


06 





47 


5 


13 


4 


94 


5 


13 


5 


35 


5 


21 


<6 


06 


<7.43 


N 7793-4 a 


359 


38553 


-32 


66690 





91 


5 


34 





46 


5 


34 


5 


32 


5 


31 


5 


67 


5 


53 


<5 


98 


<6.36 


N 7793-6 a 


359 


47568 


-32 


60091 


1 


58 


5 


14 





38 


5 


05 


4 


85 


5 


22 


5 


51 


5 


43 


<6 


38 


<7.51 


N 7793-8 a 


359 


42096 


-32 


65040 


1 


65 


5 


21 





33 


5 


20 


4 


83 


5 


19 


5 


67 


5 


24 


<6 


09 


<6.68 


N 7793-9 


359 


46133 


-32 


58006 


2 


16 


5 


05 





33 


4 


83 


4 


65 


5 


14 


5 


47 


5 


22 


<6 


94 


<7.46 


N 7793-10 


359 


41776 


-32 


61156 


2 


23 


5 


19 





32 


4 


99 


4 


68 


5 


30 


5 


61 


5 


54 


<6 


67 


<7.43 


N 7793-11 


359 


41071 


-32 


60335 


2 


48 


5 


19 





31 


4 


89 


4 


72 


5 


30 


5 


62 


5 


79 


<6 


96 


<7.50 


N 7793-12 


359 


48233 


-32 


60726 


2 


34 


5 


02 





31 


4 


73 


4 


65 


5 


06 


5 


47 


5 


54 


<6 


46 


<7.65 


N 7793-13 


359 


42245 


-32 


59289 


3 


08 


5 


23 





31 


4 


72 


4 


64 


5 


41 


5 


64 


5 


77 


<6 


93 


<7.60 


N 7793-14 


359 


47440 


-32 


57985 


2 


73 


5 


16 





30 


4 


78 


4 


61 


5 


29 


5 


59 


5 


48 


<6 


59 


<7.32 



"Identified as point sources in the 2MASS catalog l|Cutri et al.ll2003h . See Table[S]for near-IR magnitudes. 
^33-1 is Object X from lKhan et all (|201lh . 
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Table 6: Optical/NIR Luminosities of Stellar Candidates with 2MASS Counterparts 



ID 


RA 


Dec 


U 


B 


V 


R 


/ 


J 


H 


K 


M33-1 


23.35015 


30.42626 


£ 24.1 


> 24.2 


23.15 


21.61 


19.99 


17.07 


15.04 


13.6 


M33-2 


23.39209 


30.69071 


20.61 


21.84 


21.14 


20.13 


20.48 


15.96 


14.68 


14.02 


M33-8 


23.39026 


30.69038 


18.89 


19.85 


18.81 


17.94 


17.68 


16.15 


15.49 


14.25 


N 2403-1 


114.32964 


65.59473 












14.89 


14.42 


14.21 


N 2403-4 


114.22210 


65.59257 












17.21 


16.05 


14.45 


N 247-3 


11.91676 


-20.90363 












15.73 


14.79 


14.58 


N 7793-4 


359.38553 


-32.66690 












16.30 


15.59 


15.17 


N 7793-6 


359.47568 


-32.60091 












16.45 


15.98 


15.58 


N 7793-8 


359.42096 


-32.65040 












16.07 


15.56 


14.14 



IR magnitudes are presented in Table[6] along with the UBVRI magnitudes of the closest op- 
tical comrterrjarts forthe three candidates in M33 for which t here is a publicly available cata- 
log (IMassev et all 120071') . We se lect Object X (IKhan et alJboilh as a candidate (M33-1) but not 
M 33 Var A jHubble fc Sandage! Il95.4 iHumphrevs et all 1200^ )7 Although both are dust obscured 



stars with comparable bolometric luminosities and exist in our initial source list, it is apparent from 
Figure[l] that in the IRAC bands Object X is much more luminous (~ 1.5 x 10 5 L Q ) than M 33 Var A 
(~ 0.5 x 1O 5 L0). In fact M 3 3 Var A recently emerge d from a ~ 50 year dust obscured phase and 
is fading in the IRAC bands ([Humphreys et al.l 120061 ) . and we would probably have selected it in 
observations at an earlier epoch. 

It is apparent from the SEDs of the candidates (Figures[9|) that few have mid-IR luminosities 
comparable to r/Car. Moreover, the SEDs of many sources, such as N 2403-2, N 247-3 and N7793- 
10, appear so dissimilar from the SEDs of r] Car and Object X that they seem unlikely to be stellar 
sources. The MIPS band luminosity limits are also useful here. For example, the relatively low 
MIPS 70 and 160 /im luminosity limits indicate that N 247-1 must have a very flat or falling SED 
at these redder wavelengths and is probably an AGN, which is also suggested by its distance from 
the galaxy. 

On the other hand, the SEDs of M 33-5, M33-8 and M 81-10 are similar to r/Car while N300- 
1, M81-5, M81-6 and M81-7 are similar to Object X (M33-1). HST images show that M33-5, 
M33-8, and M 81-10 are compact star clusters (FigurefTOl top three panels). M81-5 is / /56 from 
a variable X-ray source with max imum luminosity of 2 x 10 38 e rgss~ 1 , which is consistent with 
the source being an X-ray binary (jRemillard McClintockl 12006). M81-7 has been classified as a 
galaxy by the SPSS su rvey but with a photometric redshift z = 0.00049 comparable to that of M81 
(jAbazaiian et al.ll2009l ) . HST images show a region of enhanced star formation consisti ng of at least 



two c omponents. Ancillary data also shows that N 7793-3 is a Hig h-Mass X- ray Binary (jMineo et al 
2012) with maximum X-ray luminosity of 3.9 x 10 37 ergss _1 ( Liul 1201 ll ). Although we do not 



currently have an explanation, 5 of the candidates in M33 (M33-3 through M33-7) are within 
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< 2'/0 of radio-selected supernova remnants ([Gordon et al.lll999l ). 



4. Rate Limits 



One advantage of searching for eruptions in the dust obscured phase is that the process is rela- 
tively easy to simulate. We eject M e of material from a star of luminosity L* and temperature T* at 
velocity v e over time period t e and assume it forms dust with total (absorption plus s cattering) visual 
opac ity Ky once it is sufficiently distant from the star. We can then use DUSTY (jElitzur fc Ivezic 



2001 



) to simulate the evolution of the mid-IR luminosities and determine the time td during which 
the source would satisfy our selection criteria. Here we use Ky = 84 cm 2 g -1 , roughly appropriate 
for silicate dust, but this is important only to the extent that the ejecta mass can be rescaled as 
M e oc Ky 1 . The key variable for estimating rates is the expansion velocity v e , be cause the detection 



period scales as td cx v e . The velocities cited for the supernova impostors (e.g., ISmith et al 



2011 



and the velocity associated with the long axis of r\ Car are high, v e > 500 kms -1 . These velocities 
are very different from those observed for the older, massive shells in the Gal axy or the shorte r 
axis of r\ Car, where v e < 100 kms -1 (see the discussion of this difference in iKochanek! 1201 lal ). 
Here we scale the results to v e = 100 kms -1 since, for example, it results in our detecting systems 
with parameters similar to rj Car at its present age, as observed, and agrees with the expansion 
velocities of the other massive Galactic shells around luminous stars. 

Detection of a shell at late times {td ^> t e ) is limited by its optical depth and temperature. 
The shell has total visual optical depth greater than ry for 

1/2 — ' M P kv \ 1/2 /100 kms" 1 



t(ry) 



M e K V 



400r, 



-1/2 



Ky 



years, 



(4) 



\-nvlry ) v \M & 100 cm 2 /g / 

and once ry < 1 it begins to rapidly fade in the mid-IR. Ignoring Planck factors, the spectral 
energy {XL\) peaks at 

he f_U N 1/4 / * « N 1/2 



A 



t 



AkT d \W 6 L Q J Vyear 100 km s" 1 J ^ ® 
so the emission peak shifts out of the IRAC bands after several decades, and our survey is primarily 
limited by the shift of the emission to longer wavelengths rather than the dec lining optical depth . 
It is better to search for these sources at 24 //m as has been done in the galaxy (Wacht er et alj|20ld : 
Gvaramadze et alj|2010i ) but that would require the resolution of JWST ([Gardner et al.l 120061 ) . A 
reasonable power-law fit to the results (—1 < log M e /M & < 1, 5.5 < log(L*/L Q ) < 6.5) of the 
DUSTY models is that the detection period is 



t d ~ t e + 66 



100 kms 



-i 



(. 



For M e ~ 10 M and L 



VlO6L 
1O 6 - 5 L like 77 Car t d g t, 



0.82 



0.043 



years. (6) 
years where t e may 

also be 50 years or more (see the discussion in lKochanek et al.ll2012al ). For present purposes, we 



190(100 kms -i/v 
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adopt td = 200 years as the period over which our selection criteria would identify an analogue of 
rj Car, consistent with the fact that our selection criteria do identify rj Car. 

We can normalize the rate of eruptions to the ccSN rate as 

(40M . \ 
-J7 N erupt RsN = fr)RsN, (7) 

M erupt J 

where Rsn is the supernova rate and all stars more massive than M erupt undergo N erupt eruptions. 



Following the rate arguments in iKochanekl (|2011al ). we can estimate the number of eruptions per 



massive star needed to explain the massive Galactic shells. If there are Ng^eii — 10 massive Galactic 
shells associated with massive stars (M > M erup t), then 

M 9 (Nsheu\ (Tv_\ V 2 / century' 1 ^ ( M erupt \ L35 / IOMq \ 1/2 / v e \ 

Irrupt - ^ 1Q J [ Qm ) ^ R SN)MW ) \ 4O M J V M e J 1 100 kms-V ' W 

where jy = 0.01 is the minimum optical depth needed to detect a shell surrounding the star and 
Rsn,mw ~ 1/century is the Galaxy's supernova rate. Since the Galactic shells are identified as 
shells primarily at 24/um, they are easier to find at low optical depths and temperatures than in 
our extragalactic survey. Thus, the massive Galactic shells imply an eruption rate relative to the 
supernova rate of > 0.2 since it is unclear whether we possess a complete inventory. Note that 
with this normalization the rate estimate does not depend on the mass scale M erup t. 



5. Conclusions 

This work empirically demonstrates that true analogs of rj Car — massive stars that have 
undergone eruptive mass ejection in the recent past (centuries) — are rare. Based on the discussion 
in SectionlU our survey can detect close analogs of r\ Car for roughly td — 200 years, consistent 
with r\ Car meeting our selection criteria. The statistics of our present sample gives us a maximum 
of N can d = 6 db 6 candidate systems after correcting for the estimated extragalactic contamination. 
Aside from the three very compact, luminous star clusters, the candidates generally do not have 
SEDs that closely resemble the SED of r^Car. Although we keep those three compact clusters in 
our candidate list for now, it is highly unlikely that they could hide luminous dusty stars similar 
to T) Car. We anticipate that further analysis and follow-up observations, using HST astrometry 
and photometry, ground based spectroscopy and Herschel 70 /xm photometry, will show that most, 
if not all, of the remaining candidates are either non-stellar or are not truly analogous to r) Car. Of 
the true stellar systems, they are clearly going to be a mixture of "eruptions" such as rj Car and 
sources with longer lived, relatively steady dusty winds such as Object X. 

We will carry out the detailed consideration of the candidates in Paper II, but suppose we scale 
conclusions about the rates to N can d = 3, which would also correspond to the 95% confidence upper 
limit we would use for estimating rate limits if we were to eliminate all the remaining candidates. 
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This implies that we are probing eruption rates of order 

d N cand / N cand \ ( 200 year ^ _ x 

-Kempt = = 0.015 I - I I — I year (9) 

for this sample of galaxies (roughly 1 per 200 years per galaxy), and fractional rates compared to 
the ccSN rate of order 

that are in the appropriate regime. In fact, it seems likely that we should not find N can d = at 
the end of Paper II, and in some senses we already have that N cand > 1 since our present sample 
contains Object X. 

Alternatively, we could estimate the expected number of candidates from ccSN rate and the 
statistics of Galactic shells as discussed in Section^ For the galaxies in our pilot study, we have 
two rather inconsistent estimates of the ccSN rate. Empirically, there were three ccSN over the 
last 20 years, which implies a rate of Rsn = 0.15 year -1 (0.05 < Rsn < 0.35 year -1 , at 90% 
confidence). On the other hand, the integrated star formation rate of the targeted galaxies implies 
a massive star formation rate, which is equivalent to the ccSN rate, of roughly Rsn = 0.014 year -1 . 
Rsn = 0.15 year -1 implies that the expected number of candidates in the targeted galaxies should 
be ~ 6 (Rsn x ft) x td for f v > 0.2 and t d = 200 years) with a > 3a chance of finding at least 1. On 
the other hand, Rsn = 0.014 year -1 reduces the probability to only about 40% and implies that 
we need to study galaxies with an integrated star formation rate of 20 Mq year -1 (10 times greater 
than what we have now) to have a > 3a chance of finding at least 1 massive dust obscured star. 

In either case, our survey can be easily expanded to at least 10 times as many galaxies (and 
integrated star formation rate) simply using archival data from the SINGS, LVL, and S 4 G surveys, 
which then probes rates far below those necessary to explain the Galactic sources. In such an 
expanded survey, additional means of suppressing contamination are important. The simplest 
method is to use the time variability of the mi d-IR emission, since ex panding shells of ejecta will 



also show a well defined pattern of fading (see iKochanek et al.l l2012al ) in the warm Spitzer bands 



(3.6 and 4.5 fj,m) and new Spitzer observations would provide a time baseline of 5-10 years to search 
for such changes. Since the principle background in our present survey appears to be extragalactic, 
time variability is a powerful means of suppressing it. Galaxies are not variable, and the mid-IR 
variability of quasars is both relatively wea k and stochastic, with a structure function of roughly 
0.1(t/4 year) 1 / 2 mag ( Kozlowski et al.ll2010i ). Two epochs separated by 6-12 months would further 



help to separate source classes by constraining variability on shorter time scales. 
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Fig. 1. — The Spectral Energy Distributions (S EP) of 77 Car now (blue solid line 



Humphreys k Davidsonlll994l). Obj ect X (black solid line, iKhan et al-lboill ). and M 33 Var A (black 
dashed line, Humphreys et al.ll2006l ). The black triangles mark luminosity at the IRAC band cen- 
ters. Although 77 Car and Object X have similar luminosities up to 3.6 /im, the SED of ryCar is 
steeply rising in the IRA C bands (a ~ 2.6; Eqn.[Tl) while Object X is almo st flat (a~ 0.2; Eqn.[T]). 



Object X, and M33 Var A ( Hubble fc Sandagelll953 ; Humphreys et al. 2006 ) are both dust obscured 



stars with comparable bolometric luminosities (jKhan et al 
is much more luminous (~ 1.5 x 1O 5 L ) than M 33 Var A 



201 lh . but in the IRAC bands, Object X 
~ 0.5 x 1O 5 L ). 
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Fig. 2. — Integrated mid-IR luminosity L m iR as a function of the slope a (Equation [lj for bright sources 
in M81. The vertical dashed lines show the slopes of blackbodies with the indicated temperatures and peak 
wavelengths (Equation EJ. The top-right (thick red) box shows the candidate selection region (L m m > 
1O 5 L and a > 0). The red triangles show the sources that also satisfy the third selection criteria, that at 
least 30% of the integrated mid-IR luminosity is emitted between 3.6 and 5.8 /xm (/ > 0.3). Of these, the 
open red triangles correspond to candidates that are known to be non-stellar in nature (see Sections 12.31 and 
13. ip . and the solid red triangles represent the surviving candidates. The green open circles show sources with 
/ < 0.3 and the black cross marks represent all the other sources. The narrow clump of points at a ~ —2.75 
correspond to normal stars with steeply falling mid-IR SEDs, while the wider clump of points to the right 
correspond to sources dominated by 8 /Ltm PAH emission. The top-left box shows the region L m iR > 10 5 ' 5 L 
and a < — 1 that was used to select normal stars in the M33 image (see Figured]). The labeled blue points 
represent objects not in M81 that are shown for comparison: Object X ("X", solid square), the compact 
cluster M33-8 ("C", open square), M33VarA (" A", large open circle), 77 Car ("77", open star), the Carina 



nebula excluding ryCar itself ("?/—", solid circle: ISmith fc Brooksl 120071 ). and the Carina nebula including 



rj Car ( u r]+" , spiked open circle; see Section HOl and Figure[6]). 
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Fig. 3. — Integrated mid-IR luminosity L m m as a function of the fraction / of L m m that is emitted 
between 3.6 and 5.8 //m for bright sources in M81. The box shows the candidate selection region 
{LmiR > 10 5 L and / > 0.3). The red triangles show the sources that also satisfy the third 
selection criteria that the mid-IR SED slope (Equations [1]) is either flat or rising (a > 0). Of these, 
the open red triangles correspond to candidates that are known to be non-stellar in nature (see 
Sections 12.31 and 13. ip . and the solid red triangles represent the surviving candidates. The green 
open circles show sources with a < and the black cross marks represent all the other sources. 
The narrow clump of points at / ~ 0.8 correspond to normal stars with steeply falling (negative 
slope) mid-IR SEDs, while the wider clump of points at / ~ 0.25 correspond to sources dominated 
by 8 jim. PAH emission. The labeled blue points are same as in FigureEJ 



Slope of SED (a) 



Fig. 4. — Extragalactic contamination for M81. Here we show all sources from a 6deg 2 region of 
the SDWFS survey transformed to the distance of M81. The symbols, lines, and axis-limits are the 
same as in Figure[2j In this SDWFS region, 449 (~ 75 deg -2 ) sources pass our selection criteria, 
indicating that we should expect ~ 13 background sources meeting our selection criteria given our 
0.17 deg 2 survey region around M81. Note that very few of the contaminating background sources 
have properties comparable to r/Car. 
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Fig. 5. — Mid and far-IR SEDs of the candidates in M33 (red lines) compared to the SEDs of 
normal stars with L m jR > 1O 55 L0, which steeply falling SEDs (mid-IR slope a < — 1, top left box 
of the Figure[2|). The dotted portions of the SEDs correspond to the MIPS 70 and 160 ^m flux 
upper limits. The SED of Object X is highlighted (red-black lighter dashed line) and r/Car (black 
heavier dashed line) is shown for comparison. 
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Fig. 6. — The SEDs of rj Car ( "ry" , black tria ngles , Humph reys Davidsonlll994l ). the Carina nebula 
excluding rj Car itself ( "ry-" , blue squares, ISmith &: Brooks! 120071 . and the entire dusty complex 
containing rj Car and other massive stars including rj Car ("77+", red circles, Section [2^|) . The first 
two SEDs are spline interpolated and summed to produce the third. The SED of the compact 
cluster M33-8 ("C", green dashed line, HST image in FigurefTUj) is shown for comparison. In 
Figures El El and H] we label these rj, 77—, rj+, and "C" respectively. 
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Fig. 7. — SEDs of four different classes of objects that met our selection criteria: a candidate dusty 
star in NGC2403, a star-cluster in M33, a QSO behind M81, and a galaxy behind NGC 7793. 
FiguredU] shows IRAC and HST images of the compact cluster and the galaxy. 
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Fig. 8. — SEDs of sources that met our selection criteria but were rejected due to association with 
non-stellar sources. The dotted portions of the SEDs correspond to the MIPS 70 and 160 /im flux 
upper limits. The SED of r] Car (dashed blue line) is shown for comparison. 
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Fig. 9. — SEDs of sources that met our selection criteria and were not rejected due to association 
with non-stellar sources. The dotted portions of the SEDs correspond to the MIPS 70 and 160 fim 
flux upper limits. The SEDs of rj Car (dashed blue line) and Object X (dot-dashed black line) are 
shown for comparison. 




Fig. 10.— IRAC and HST images of the compact stellar clusters M33-5, M33-8, and M 81-10, 
and the background galaxy N 7793-2. The clusters are resolved in the HST images with FWHM of 
/ /87 ~ 4.1 pc (M 33-5), 0"77 ~ 3.6 pc (M 33-8) and 0"34 ~ 6.1 pc (M 33-8). They are very luminous 
(fewxlO 7 L0) and their SED shapes are very similar to ryCar (Figure[9|. 



